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Abstract 

We study the axial amplitudes for the N-A transition in models with quarks and 
chiral mesons. A set of constraints on the pion field is imposed which enforces PCAC 
and the off-diagonal Goldberger-Treiman relation. The quark contribution to the 
amplitudes in general strongly underestimates the amplitude as well as the vrNA 
strong coupling constant. We show that the results are considerably improved in 
models that, in addition to the pion cloud, incorporate a fluctuating a field inside 
the baryon. 
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1 Introduction 



The structure of the weak axial N-A transition currents is ideally probed in 
neutrino or charged- lept on scattering experiments on deuterium or hydrogen. 
The experimental efforts so far have been focused on the determination of the 
dipole axial mass parameter [1], without an attempt to break down the tran- 
sition current into form- factors [2]. Although a number of phenomenological 
predictions for the dominant coupling C^(0) exist (see Table I of [3] for an 
exhaustive list), the dependence of the form-factors on momentum transfer is 
very poorly known. Data on the non- leading form-factors C^{Q^) and C^{Q'^) 
are especially scarce [4] . New information on the weak axial form-factors is ex- 
pected from parity-violating electron scattering experiments planned at Jef- 
ferson Laboratory [5]. 

Theoretical investigation of axial transition amplitudes in different versions of 
the quark model is of particular interest since it may reveal the importance of 
non-quark degrees of freedom in baryons, in particular the chiral mesons. Yet, 
except for the calculation in the non-relativistic quark model [6], there exist 
almost no model predictions for the axial transition amplitudes. This can be 
traced back to the difficulty of incorporating consistently the pion field which 
is necessary to describe the proper low-Q^ behaviour of the amplitudes. 

The lack of experimental and theoretical knowledge in the weak sector is 
in contrast to the case of electro-excitation of the A resonance, which has 
been extensively studied theoretically in the constituent quark models [7] as 
well as chiral models [8], and experimentally [9]. In [10] we have pointed out 
the important role played by the pion cloud in the determination of electro- 
production amplitudes, in particular to the E2/M1 and C2/M1 ratios. This 
has later been confirmed in other chiral models [11] and dynamical approaches 
[12]. 

The aim of this work is to study some general properties of the axial ampli- 
tudes in chiral quark models and present theoretical predictions in two typical 
representatives of such models, the linear a model with quarks and the Cloudy 
Bag Model. We derive a set of constraints on the pion field which enforce the 
proper behaviour of the amplitudes in the vicinity of the pion pole. We also 
address the long-standing problem of a too low ttNA coupling constant which 
rather systematically appears in all quark models. Comparing the results in 
the two models we are able to draw some general conclusions regarding the 
contribution of chiral mesons to the weak amplitudes as well as to the strong 
ttNA form factor. 
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2 The cixial transition amplitude and the ofF-diagonal Goldberger- 
Treiman relation 



The axial N-A transition ampUtude is usually parameterized in terms of the 
Adler form-factors C^{Q'^) as ^ 



where p'^ = (Ma; 0,0,0), is the corresponding Rarita-Schwinger spinor, 
p is the four-momentum of the nucleon and = (A'o; 0, 0, k) is the four- 
momentum of the incident weak boson. Then — k'^ = = —Q^ and k^ — 
{Ml - - Q2)/2Ma. For simplicity, we take the third isospin component 
(a = 0) of the axial current. We have omitted from (1) the C^{Q^) term 
[2], which is consistent with the prediction of quark models in which quarks 
occupy only the I — state. 

It is convenient to work with helicity amplitudes^ 



^^ = -(A+(pO,«a = i I A',iO) I pb),SN = i) , (2) 

if = -(A+(p'), «A = 1 1 • A%0) I p(p), = I) , (3) 

if = -(A+(p'), «A = i I • A0(0) I Pip), = -i) , (4) 

- -(A+(p'), SA = II £o ■ A°(0) I p(p), SN = I) , (5) 

where s denotes the third spin component, and e are the usual polarisation 
vectors. The helicity amplitudes are related to the Cf^ form-factors by 



c 



Ml 




/3 (ko^A ^0 f A 



+ 



/r2 



k^ . 



(6) 
(7) 



^ Definition of the transition current with respect to the brings an additional 
isospm factor V3 to RHS of (1). 

^ The hehcity amplitudes arc normally defined as the matrix elements of the in- 
teraction Hamiltonian and contain an additional factor y/4:7raw f^Ko, e. g. = 
\J Airaw 1 2i^o ■, where Kq = ko{Q^ = 0) and aw is the weak fine-structure con- 
stant. 



3 



Ml 




(8) 



In the approximation with = we have only one independent transverse 
amphtude, since in this case ^4^ = y/SAj^. 

From (1) it follows that the divergence of the transition axial current vanishes 
in the chiral limit provided Cq{Q'^) = MlC^{Q'^)/Q'^. The pole behaviour of 
the Cq amplitude suggests that it is related to the term in the axial current 
responsible for the pion decay, ^a(poie)(^) ~ f-K&^T^aix), where /^r = 93MeV is 
the pion decay constant. Taking a finite mass for the pion the divergence does 
not vanish but is replaced by PC AC: 



(A+(pO I I N+(p)) = -ml A(A+(p') 1 7r,(0) | N+(p)) , 



(9) 



where the transition matrix element of the pion field is related to the strong 
form factor G'^an(Q^) by 

(A^M|..(0)|Nnrt)^i^^^/|. (10) 

Assuming that ^a(poic)(^) dominates the amplitude for — > —ml, we 
obtain the off-diagonal Goldberg er-Treiman relation [2,13,14]: 

Q^^-^l. (11) 

For a smooth interpolating pion field wc expect that (11) holds also for mod- 
erate Q"^ in the physically accessible region. 



3 Helicity amplitudes in chiral quark models 



For a variety of models involving quarks interacting with chiral fields a and tt 
the Hamiltonian can be written in the form 



H 



+ (V^ + mDn^] + U{a, n) + Y.ja'Ka ,(12) 



where ja is the quark source, P^^ is the pion conjugate momentum, and 
H„ are the free-quark and the cr-meson terms, and U (cr, tt) is the meson self- 
interaction term. In the Cloudy Bag Model the a field and the self-interaction 
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term arc absent, while in the hnear a-modcl all terms are present and the 
self- interaction term is the Mexican-hat potential (see (22) below). 

Let |N) and | A) be the exact solution of the Hamiltonian for the ground state 
and for the A, respectively, with H\N) ^ En\N) and //|A) ^ Ea\A). Then 
(N|[i/,P^]|N) = (A|[F,P^]|A) = and (A| [if, P^] |N) = i(PA - ^N)'(A|7f|N). 
Evaluating the commutators using (12) for a = 0, we obtain 



(-A + ml){N 1 7ro(r) | N) = -(N | Mr) \ N) , (13) 
(-A + ml) {A 1 7ro(r) | A) = - (A | Mr) \ A) , (14) 
(-A + ml- (Ea - E^)'){A 1 7ro(r) | N) = - (A | Mr) I N) . (15) 

The sources on the RHS of (13)-(15) consist of the quark term and the term 
originating from the meson self-interaction (if present) : 



These relations hold for the exact solutions of (12). In an approximate com- 
putational scheme they can be used as constraints. 

We now show an important property of the axial transition amplitudes be- 
tween states which satisfy these virial relations. We split the axial current into 
the non-pole and the pole part, A" — ^°non-poie) + ^(poie)' where 



^fnon-pole) = V^T^TS '^T^ + {t ' A)^"^ " , (17) 

^(pole)=A5"^- (18) 

Since the pole part involves only the pion field we can use (15) to evaluate its 
contribution to the amplitudes. Note that (15) is equivalent to (10) since in 
our model we can write the strong N-A transition form-factor as 

^f^4(A||.o(0)||N). (19) 



We find if (pole) = and 



qA _kofA _ 2 G^NAjQ^) Ukko 

^ipoie) - ^ ^(poie) - 3 Q^ + ml' ^^^^ 
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The pole term (18) contributes only to Cq, 



iA _ f <^7rNA_(Q 




<^6(pole) - /tt ^2 , _2 \/ Q ' y"^^) 



while C*4(poie) = ^5\poie) ~ ^- conclude that in models in which the pion 
contribution to the axial current has the simple form f.^d'^'Ka and the pion 
field satisfies the virial relation (15) there is no pion contribution to the 
and amplitudes while Cq is almost entirely dominated by the pion pole. 
In such models only the quarks contribute to the and amplitudes. In 
this respect, the calculation of in a constituent quark model calculation 
(e. g. [6]), is still legitimate. 



4 Constrained calculation in the linear cr-model 

The hnear cr-model assumes the following form of jt and U [15]: 

it = i 5 E qa^mi , U^j{a' + 7r'- flf . (22) 

i=l ^ 



Here is the quark bispinor for the valence orbit (assumed to be different 
for the nucleon and the A), and A = (m^ — m^)/2/J. The free parameters 
of the model are the couphng strength g related to the "constituent" mass 
of the quark gf^^, and the mass of the a meson m^. The model has been 
successfully applied to the description of the nucleon and A properties. So far 
the physical states have been constructed from the mean-field solution using 
either cranking [16] or the Peierls-Yoccoz projection [17]. In the latter method 
the mean-field solution for the pion field is interpreted as a coherent state. 

The mean-field solution fulfills the diagonal virial relations (13)-(14) but not 
the off-diagonal relation (15). To satisfy this relations it is necessary to include 
a channel representing the A decay, i.e. a term that asymptotically represents 
the nucleon and a free pion. We have therefore taken a more general ansatz 
for the A: 

|A) =A/a{p^$a|A,) + J dM(fc)[«L(fc)|N)]^^} , (23) 



where the first term represents the bare A state surrounded by a cloud of pions 
and a mesons, is the projection operator on the subspacc with isospin and 
angular momentum |, |N) is the nucleon ground state, and [ ]^^ denotes a 
pion-nucleon state with isospin | and spin |. Requiring that the energy of this 
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state is stationary, the denominator of ri{k) takes the form ouk — {E^ — E-^) 
which is also the form imphed by (15). For the nucleon we assume: 

|N) = Af^P'^ [$n|N,) + $na| A,)] . (24) 

Here <I>n and $na stand for hedgehog coherent states describing the pion 
cloud around the bare nucleon and bare A, respectively. To match the third 
constraint, (15), the denominator of the pion state in the second term of (24) 
should behave as ujk + wq with ujq = [E^ — -En). In the above ansatz, only one 
profile for the a field is assumed ^ . 

The properties of the ground state are dominated by the first term in (24), and 
imposing the off-diagonal constraint influences only slightly the results. For 
the A, the inclusion of the decaying channel modifies the long-range behaviour 
of the pion field, and yields the correct low-Q^ behaviour of the transition 
amplitudes as explained in the previous section. The calculated A-N splitting 
is typically only (50 — 70) % of the experimental value. In order to make a 
sensible comparison of the transition amplitudes with the experimental ones, 
it is necessary to have the correct kinematical relations in the model. This 
can be achieved by including an additional phenomenological term in the 
Hamiltonian mimicking either the chromo-magnetic or the instanton-induced 
interaction between quarks and adjusting its strength such as to bring the 
A-N splitting to the experimental value. 



5 Calculation of the amplitudes 

We calculated the amplitudes in two models: in the linear cr-model and in the 
Cloudy Bag Model. In the Cloudy Bag Model we assume the usual pertur- 
bative form for the pion profiles [18] using the experimental masses for the 
nucleon and A, which fulfills the virial constraints (13)-(15). Since the pion 
contribution to the axial current in the Cloudy Bag Model has the form of the 
pole term in (18), only the quarks contribute to the and amplitudes. 

The amplitudes (2)- (5) are defined between states with good 4-momenta p' 
and p respectively while in the model calculations localised states are used. 
We can use such states in our calculation of amplitudes by interpreting them 
as wave packets of states with good linear momentum. Extending the method 
explained in [14] we find for a chosen component of the axial current evaluated 
between locahsed states, (A|74(r)|N): 

I dpipl{p + k)ip^{p){A{p + k)\A{0)mp)) = I dre'^^iAlAirm .(25) 
^ Since the a field is scalar its analog of (15) is identically zero. 
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Here the matrix element of (1) is taken on the LHS and '■Pn{p) and v'a(p) are 
(normahsed) functions describing the center-of-mass motion of the locahsed 
solution for the nucleon and the A, respectively. We assume that the spread 
of the wave packet is of the order of the inverse baryon mass (M~^) and 
use for simplicity the same spread for the nucleon and the delta. The Adler 
form-factors of (6)-(8) are then modified in such a way that and are 
multiplied by the factor 2Ma/(Ma + Mn), while 



kok Ma + Mis 




Ml 

2M^ ^ 



(26) 



We have neglected terms of the order fc^/M^. Similarly, the strong Cttna form 
factor (19) acquires the same correction factor. The essential property that 
the pole contribution cancels out in and still persists as well as the 
relation (21) for Cq. 




Q- [ (GeV/c)^ 



Fig. 1. The amplitude C^{Q'^) in the linear ir-model. The experimental value of 
1.22 ± 0.06 at Q2 = [21] is based on data from ANL and BNL [22,23]. The 
error ranges are given by the spread in the axial-mass parameter M\ as determined 
from neutrino scattering experiments (broader range) and from electro-production 
of pions [1] (Ma = (1.077 ± 0.039) GeV, narrower range). Full curves: wave-packet 
result; dashed curves: calculation from Gt^na (H)- 

Fig. 1 shows the amplitude in the linear cr-model with g = 4.3 and 
m„ = 600 MeV compared to the experimental weak axial form-factors given 
in the convention of Adler [19,20], with a phenomenological dipole parameter- 
isation C^iQ^) = C^(0)/(1 + Q'^/Mlf. The ^^(0) is 25% higher than the 
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experimental estimate, while the Ma from a dipole fit to our calculated values 
matches the experimental Ma to within a few percent. 

We note that for the nucleon we obtain = 1-41 which is roughly the same 
amount higher than the experimental value of 1.27. On the other hand, if 
we determine C^{Q^) from the calculated strong ttNA form-factor using the 
Goldberger-Treiman relation (11) we obtain a better agreement, yet with a 
steeper fall-off corresponding to Ma ~ 0.80 GeV. The discrepancy between the 
two calculated values (17% at Q'^ — — where (11) holds) is a measure for 
the quality of our approximate computational approach. It can be attributed 
to a too large meson contribution originating from the last two terms in (17). 
Since in this model only the meson fields bind the quarks it is reasonable 
that their strength is overestimated in the variational calculation. The effect 
of the meson self-interaction (the second term in (16)) is relatively less pro- 
nounced in the strong coupling constant (only ~ 20%) than in C^{Q'^). Both 
G7rNA(0) and G,rNN(0) arc over-estimated in the model by ~ 10%. Still, the 
ratio G,rNA(0)/G7rNN(0) = 2.01 is considerably higher than either the familiar 
SU(6) prediction ^72/25 or the mass-corrected value of 1.65 [14], and com- 
pares reasonably well with the experimental value of 2.2. This improvement 
is mostly a consequence of the renormalisation of the strong vertices due to 
pions. 

The value of grows with g and rria in contrast to Gttna which remains 
almost constant over a large range of model parameters. In our calculation 
we cannot use much lower values for g since the solution becomes numerically 
unstable. 

In the Cloudy Bag Model the picture is reversed. Here only the first term 
in (17) contributes to the amplitudes; as a result the amplitude is less 
than 2/3 of the experimental value (see Fig. 2). The behaviour of is 
similar as in the pure MIT Bag Model (to within 10%), with fitted Ma ~ 
1.2GcVfm/i?. The off-diagonal Goldberger-Treiman relation is satisfied in 
the Gloudy Bag Model, but from Gt^na has a steeper fall-off with fitted 
Ma ~ 0.8GeVfm/i?. The ratio C^{0)/gA is close to the model-independent 
prediction of [24]. 

The large discrepancy can be to some extent attributed to the fact that the 
Cloudy Bag Model predicts a too low value for GttNN, and consequently GttNA- 
Taking a smaller value of /tt in order to increase the strong coupling constants 
does not improve the results since on the RHS of (11) compensates for the 
change in GttNA- We have found that the pions increase the Gttna/Gttnn ratio 
by ~ 15 % through vertex renormalisation. The effect is further enhanced by 
the mass-correction factor 2Ma/(Ma-|-Mn), yet suppressed in the kinematical 
extrapolation of G7rNA(Q^) to the SU(6) limit. This suppression is weaker at 
small bag radn R: the ratio drops from 2.05 a,t R — 0.7 fm to 1.60 (below the 
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Fig. 2. The amplitude C^{Q'^) in the Cloudy Bag Model for three values of the bag 
radius. Experimental uncertainties are as in caption to Fig. 1. 

SU(6) value) at R = 1.3 fm. 

The determination of the is less reliable because the meson contribution 
to the scalar amplitude is very sensitive to small variations of the profiles. 
However, the experimental value is very uncertain as well. Neglecting the 
non-pole contribution to 5*^ and Cq (the pole contribution cancels out) we 
see from (26) that the value of is dominated by the term —{M^/2M'^) C^, 
in agreement with the popular parameterisation of the amplitudes. In our 
models, the non-pole contribution to Cq is not negligible and tends to increase 
at small Q^, as seen in Fig. 3. 

The amplitude is governed by the pion pole for small values of and 
hence by the value of GttNA which is well reproduced in the linear a-model, 
and underestimated by ~ 35 % in the Cloudy Bag Model. Fig. 4 shows that 
the non-pole contribution becomes relatively more important at larger values 
of Q2. 



6 Conclusions 

To the best of our knowledge the present work is the first attempt to calculate 
the axial N-A transition amplitudes in a quark model which consistently in- 
cludes the chiral mesons already at the Lagrangian level. We have derived a set 
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Fig. 3. The amplitude C^{Q'^) in the hnear cr- model, with model parameters and 
experimental uncertainties due to the spread in Ma as in Fig. 1. Experimentally, 
^4(0) = —0.3 lb 0.5 [4]. For orientation, the value for C^{0) is used without er- 
ror-bars. 

of constraints which ensures the proper treatment of the pion pole dominating 
the transition at low Q^. Though there is a rather strong discrepancy between 
calculated amplitudes in the two models considered here, we are nonetheless 
able to draw some general conclusions about the role of the chiral mesons. 

The quark contribution alone strongly underestimates the amplitude. 
Models in which only a linear coupling of pions to quarks is added do not 
improve the situation since in such the pion term in the axial current 

does not contribute to the amplitude. On the other hand, the inclusion of 
meson self-interaction which allows for a substantial deviation of the a field 
from its vacuum value inside the baryon considerably increases C^. The lin- 
ear cr-model seemingly overestimates this contribution as it could have been 
anticipated from the overestimate of qa obtained in this model. 

Regarding the ratio GTrNA/G^rNN we find that it is the pion cloud which en- 



hances its value compared to the SU(6) value of y 72/25; in the linear cx-model 
as well as in the Cloudy Bag Model for smaller bag radii the ratio is greater 
than 2 and not far from the experimentally determined value of 2.2. 

The Q^-behaviour of the axial amplitudes is well reproduced in the linear cr- 
model. We stress that the behaviour of G,rNA(Q^) is considerably softer, with 
a cut-off parameter (corresponding to the axial mass Ma) of ~ 0.8 GeV. A 
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Fig. 4. The non-pole part and the total amplitude C^(Q^) in the linear cr-model. 

Model parameters are as in Fig. 1. 

similar trend is also seen in the Cloudy Bag Model for bag radii above ~ 1 fm. 
The popular assumption in which the same value for the strong and axial 
cut-offs are taken is therefore not supported by the two models. 

This work was supported by FCT (POCTI/FEDER), Lisbon, and by The 
Ministry of Science and Education of Slovenia. 
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